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Minimal Actuator Placement with Optimal Control Constraints 


V. Tzoumas, M. A. Rahimian, G. J. Pappas, A. Jadbabaie* 


Abstract —We introduce the problem of minimal actuator 
placement in a linear control system so that a bound on the 
minimum control effort for a given state transfer is satisfied 
while controllability is ensured. We first show that this is an 
NP-hard problem following the recent work of Olshevsky QJ. 
Next, we prove that this problem has a supermodular structure. 
Afterwards, we provide an efficient algorithm that approximates 
up to a multiplicative factor of O(logn), where n is the size of 
the multi-agent network, any optimal actuator set that meets the 
specified energy criterion. Moreover, we show that this is the best 
approximation factor one can achieve in polynomial-time for the 
worst case. Finally, we test this algorithm over large Erdos-Renyl 
random networks to further demonstrate its efficiency. 

Index Terms —Multi-agent Networked Systems, Transfer to a 
Single State, Leader Selection, Minimal Network Controllability. 

I. Introduction 

During the past decade, control scientists have developed 
various tools for the regulation of large-scale systems, with 
the notable examples of El for the control of biological sys¬ 
tems, El for the regulation of brain and neural networks, 01 
for network protection against spreading processes, and El 
for load management in smart grid. On the other hand, the 
enormous size of these systems and the need for cost-effective 
control make the identification of a small fraction of their 
nodes to steer them around the state space a central problem 
within the control community iH, [|6l, El. El. (SI. ifTOl . 

This is a combinatorial task of formidable complexity; as 
it is shown in El. identifying a small set of actuator nodes 
so that the resultant system is controllable alone is NP- 
hard. Nonetheless, a controllable system may be practically 
uncontrollable if the required input energy for the desired state 
transfers is forbidding, as when the controllability matrix is 
close to singularity ED- Therefore, by choosing input nodes 
to ensure controllability alone, one may not achieve a cost- 
effective control for the involved state transfers. In this paper, 
we aim to address this important requirement, by introducing 
a best-approximation polynomial-time algorithm to actuate 
a small fraction of a system’s states so that controllability 
is ensured and a specified control energy performance is 
guaranteed. 

In particular, we consider the selection of a minimal number 
of actuators such that a bound on the minimum control effort 
for a given transfer is satisfied while controllability is ensured. 
Finding the appropriate choice of such a subset of nodes is a 
challenging task, since the search for a subset satisfying certain 
criteria constitutes a combinatorial optimization problem that 
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can be computationally intensive. Indeed, it is shown in El 
that identifying the minimum number of actuators for inducing 
controllability alone is NP-hard. Therefore, we extend this 
computationally hard problem by imposing an energy con¬ 
straint on the choice of the actuator set, and we solve it with 
an efficient approximation algorithm. 

Specifically, we first generalize the involved energy ob¬ 
jective to an e-close one, which remains well-defined even 
when the controllability matrix is non-invertible. Then, we 
make use of this metric and we relax the controllability 
constraint of the original problem. Notwithstanding, we show 
that for certain values of e all solutions of this auxiliary 
program still render the system controllable. This fact, along 
with a supermodularity property of the generalized objective 
that we establish, leads to a polynomial-time algorithm that 
approximates up to a multiplicative factor of O(logn) any 
optimal actuator set that meets the specified energy bound, 
when the latter lies in a certain range with respect to n. 
Moreover, we show that this is the best approximation factor 
one can achieve in polynomial-time for the worst case. 

Hence, with this algorithm we aim to address the open 
problem of actuator placement with energy performance guar¬ 
antees HI, El, El, EH, ini. To the best of our knowledge, 
we are the first to study the selection of a minimal number 
of actuators so that a bound on the minimum control effort 
for a given transfer is satisfied. Our results are also applicable 
to the case of average control energy metrics El and can be 
extended to the cardinality-constrained actuator placement for 
minimum control effort, where the optimal actuator set is se¬ 
lected so that these metrics are minimized, while its cardinality 
is upper bounded by a given value. These and other relevant 
extensions are explored in the companion manuscript EOl. 

The remainder of this paper is organized as follows. The 
formulation and model for the actuator selection problem are 
set forth in Section m In Section HII] we discuss our main 
results, including the intractability of this problem, as well as 
the supermodularity of the involved control energy objective. 
Then, we provide an efficient approximation algorithm for its 
solution. Finally, in Section |IV] we illustrate our analytical 
findings on an integrator chain network and we test their per¬ 
formance over large Erdos-Renyi random networks. Section fV] 
concludes the paper. All proofs can be found in the Appendix. 

II. Problem Formulation 

Notation: Denote the set of natural numbers {1, 2,...} as 
N, the set of real numbers as K., and let [n] = {1, 2,..., n} for 
all n € N. Also, given a set X, denote \X\ as its cardinality. 
Matrices are represented by capital letters and vectors by 
lower-case letters. For a matrix A, is its transpose, and 
Aij is its element located at the *—th row and j—th column. 


For a symmetric matrix A, A = A"^; and if A is positive 
semi-definite, or positive definite, we write A ^ 0 and A 0, 
respectively. Moreover, for i € [n], let be an n x n matrix 
with a single non-zero element: la = 1, while Ijk = 0, for j, 
k ^ i, and denote the identity matrix by I, where its dimension 
is inferred from the context. Additionally, for 6 G K", diag{5) 
denotes an n x n diagonal matrix such that diag{S)ii = Si, 
for all i G [n]. 


A. Actuator Placement Model 

Consider a linear system of n states, xi,X 2 , ■ ■ ■, Xn, whose 
evolution is described by 

x{t) = Ax{t) + Bu{t),t > to, (1) 

where fo G R is fixed, x = {x\,X 2 , ■ ■ ■ ,Xn}, x{t) = dx/dt, 
while u is the corresponding input vector. The matrices A and 
B are of appropriate dimension. Without loss of generality, 
we also refer to O as a network of n agents, 1, 2,..., n, 
which we associate with the states Xi,X 2 t ■ ■, Xn, respectively. 
Moreover, we denote their collection as V = \n]. Henceforth, 
the interaction matrix A is fixed, while a special structure is 
assumed for the input matrix B. 

Assumption B = diaglS), where S G {0,1}". 

Each choice of the binary vector 5 in Assumption[T]signifies 
a particular selection of agents as actuators. Hence, if Si = 1, 
state i may receive an input, while if Si = 0, receives none. 
We collect the above and others into the next definition. 

Definition 1 (Actuator Set, Actuator). Given S G {0,1}" and 
B = diag{S), let IS. GV be such that \/i G A, Si = 1, while 
Vi ^ A, Si = 0; then, A is called an actuator set and any 
agent i G A is called an actuator. 


B. Controllability and the Minimum Energy Transfer Problem 

We consider the notion of controllability and relate it to the 
problem of selecting a minimum number of actuators for the 
satisfaction of a control energy constraint. 

Recall that ([U is controllable if for any finite ti > to and 
any initial state xq = x{to), the system can be steered to any 
other state xi = x{ti), by some input u{t) defined over [fo, fi]. 
Moreover, for general matrices A and B, the controllability 
condition is equivalent to the matrix 

r{to, h)= [" dt, ( 2 ) 

Jto 

being positive definite for any ti > to HI]. Therefore, we 
refer to r(fo,fi) as the controllability matrix of ([U. 

The controllability of a linear system is of great interest, 
because it is related to the solution of the following minimum 
energy transfer problem 


/•ti 

minimize / u{t)'^u{t) dt 

“(•) Jto 

subject to (3) 

x{t) = Ax{f) + Bu{t), to <t<ti, 
x(to) = xo,x(ti) = Xi, 

where A and B are any matrices of appropriate dimension. 
In particular, if ([T]i is controllable for the given A and B, the 
resulting minimum control energy is given by 

{xi - e'^^a;o)^r(fo, ti)~^{xi - e^'^xo), (4) 

where t = ti — to lfT4l . Therefore, if xi — e^'^xo is spanned 
by the eigenvectors of T{to,ti) corresponding to its smallest 
eigenvalues, the minimum control effort (|4]i may be forbid¬ 
dingly high HD. Hence, when we choose the actuators of a 
network so that controllability is ensured and an input energy 
constraint for a specified state transfer is satisfied, we should 
take into account their effect on 

Moreover, controllability is an indispensable property for 
any linear system, while in many cases is viewed as a structural 
attribute of the involved system na that holds true even by 
any single input nodes, as in large-scale neural networks 0. 
This motivates further the setting of this paper, where the 
actuators are chosen so that a bound on the minimum control 
effort for a given transfer is satisfied and overall controllability 
is respected. 

Per Assumption HI some further properties for the controlla¬ 
bility matrix are due. First, given an actuator set A, associated 
with some S, let Fa = r(fo,fi); then, 

n 

rA = ^<5,r„ (5) 

where for any i G [n], F^ = * dt, that is, each 

Fi is a constant positive semi-definite matrix determined by 
A, to and ti. To see why (l5]l holds true, observe that B = 
diag{S) implies B — BB^ = and Q follows 

upon replacing this in (l2|l. Furthermore, note that (l5]l together 
with the fact that F^ ^ 0, for any i G [n] gives Fai A Faj 
whenever Ai C A 2 . 

C. Actuator Placement Problem 

We consider the problem of actuating a small number of 
system’s o states so that the minimum control energy for 
a given transfer meets some specified criterion and controlla¬ 
bility is ensured. The challenge is in doing so using as few 
actuators as possible. This is an important improvement over 
the existing literature where the goal of actuator placement 
problems have either been to ensure just controllability HI 
or the weaker property of structural controllability na, HD- 
Other relevant results consider the task of leader-selection na, 
HD, where the leaders, i.e. actuated agents, are chosen so as 
to minimize an appropriate mean-square convergence error of 
the remaining agents. Our work also departs from a set of 
works that study average energy metrics, such as the minimum 


eigenvalue of the controllability Gramian or the trace of its 
inverse 0, a, 113. Instead, here we consider an exact energy 
objective and require it to satisfy a particular upper bound. 

Let Cr = {A C V : |A| < r, La 0} be the actuator sets 
of cardinality at most r that render ([Til controllable. Then, for 
any A C V, we write A G C|a| to denote that A achieves 
controllability. Furthermore, we set 

V = {xi - e^'"xo)/\\xi - e^'^xo\\2. 

We consider the problem 

minimize IAI 
ACV 

subject to 

A G C|A|) 

v^T'^^v < E, 

for some positive constant E. This problem is a generalized 
version of the minimal controllability problem considered 
in m, so that its solution not only ensures controllability, 
but also provides a guarantee in terms of the minimum input 
energy required for the normalized transfer from xq to xi\ 
indeed, for E ^ oo, we recover the problem of Cl. 

For some extra properties of dD, note that for any A G C|a|, 
0 -< Fa ^ Fy, i.e. ll20l . Hence, ([III is 

feasible for any E such that 

< E. (6) 

Observe that this lower bound depends only on A and v, i.e. 
also on n, as well as on tg and ti. 

Moreover, dD is NP-hard, since it looks for a minimal 
solution and so it asks if C^. ^ 0 for any r < n m. Thus, 
we need to identify an efficient approximation algorithm for 
its solution, which is the subject of the next section. 

III. Minimal Actuator Sets with Constrained 
Minimum Energy Performance 

We present an efficient polynomial-time approximation al¬ 
gorithm for dD- To this end, we first generalize the involved 
energy objective to an e-close one, that remains well-defined 
even when the controllability matrix is non-invertible. Next, 
we relax (|D by introducing a program that makes use of 
this objective and ignores controllability constraint of (|D. 
Nonetheless, we show that for certain values of e all solutions 
of this auxiliary program still render the system controllable. 
This fact, along with the supermodularity property of the 
generalized objective that we establish, leads to our proposed 
approximation algorithm. The discussion of its efficiency ends 
the analysis of dD- 

A. An e-close Auxiliary Problem 

Consider the following approximation to Problem (|D 

minimize IAI 
ACV 

subject to ) 

m) < E, 


where (j){A) = v'^{VA + el) + v'[{rA + e'^I) ^Vi, 

for any ACV, while vi,V 2 ,. ■. ,Vn-i are an orthonormal 
basis for the null space of v, and e is fixed such that 0 < e < 
1/E, given E. Observe that the controllability constraint is 
now ignored, while the energy objective is well-defined for any 
actuator set A including the empty set, since the invertibility 
of Fa -f e/ and Fa -f e^I is always guaranteed for e > 0. 

The e-closeness is evident, since for any A G C|Ap </>(A)^ 
as e —>■ 0. Notice that we can take e —>• 0, since we 
assume any positive e < 1/E. 

B. Approximation Algorithm for Problem dB 

We first prove that all solutions of dD for 0 < e < 1/E, 
render the system controllable, notwithstanding that no con¬ 
trollability constraint is imposed by this program on the choice 
of the actuator sets. Moreover, we show that the involved e- 
close energy objective is supermodular, and then we present 
our approximation algorithm, followed by a discussion of its 
efficiency, which ends this subsection. 

Proposition 1. Fix cv > 0. Then, Ve, 0 < e < l/w, ;/VA C V, 
4>[A) < Lo, then A G C|a|. 

Note that uj is chosen independently of the parameters 
of system dD- Therefore, the absence of the controllability 
constraint at Problem dD for 0 < e < l/i?is fictitious; 
nonetheless, it obviates the necessity of considering only those 
actuator sets that render the system controllable. 

The next lemma is also essential and suggest an efficient 
approximation algorithm for solving dD- 

Proposition 2 (Supermodularity). The function v'^{Ta + 
el)~^v -h vj {1^ A + e^I)~^Vi : A C V i-A R A 

supermodular with respect to the choice of A. 

Inspired by the literature on set-covering problems subject to 
submodular constraints, ED, El, El, we have the following 
efficient approximation algorithm for Problem dD, and, as we 
illustrate by the end of this section, for Problem (jD as well. We 
note that a corollary of the above proposition is that is 

supermodular as well, but over the sets A C V that render dD 
controllable. 


Algorithm 1 Approximation Algorithm for the Problem (jp. 
Input: Upper bound E, approximation parameter e < 1/E, 
matrices Fi, r2,..., r„, vector v. 

Output: Actuator set A. 

A4-0 

while (/'(A) > E do 

Qi G argmax^gv\A{'i^(^) “ U {a})} 

A i — A U 

end while 


For the efficiency of Algorithm [D the following is true. 

Theorem 1 (A Submodular Set Coverage Optimization). 
Denote as I* the cardinality of a solution to Problem (jD <^nd 





as A the selected set by Algonthm\I\ Then, 


^ S C\A\, 

(/.(A) < E, 

|A| 

F = 0(log n + log €~^ + log 


ne ^ — <?!>(V) 


= F, 


1 


E - ^(V) 


)• 


(7) 

( 8 ) 

(9) 

( 10 ) 


Therefore, the polynomial-time Algorithm [T] returns a set of 
actuators that meets the corresponding control energy bound 
of Problem dFl i. while it renders system ([T]i controllable. More¬ 
over, the cardinality of this set is up to a multiplicative factor 
of F from the minimum cardinality actuator sets that meet 
the same control energy bound. In Section IIII-CI we elaborate 
further on the dependence of this multiplicative factor on 
n, e and E, using (fTOl l. while in Section IIII-DI we finalize 
our treatment of Problem (HJl by employing Algorithm [T] to 
approximate its solutions. 


C. Quality of Approximation of Algorithm\I\for Problem (|F]l 

The result in ( fTOb was expected from a design perspective: 
Increasing the network size n or improving the accuracy by 
decreasing e, as well as demanding a better energy guarantee 
by decreasing E, should all push the cardinality of the selected 
actuator set upwards. Also, note that loge“^ is the design 
cost for circumventing the difficulty to satisfy controllability 
constraint of Problem (U) directly HI . 

Furthermore, per (doll and with E — ^(V) and e both 
fixed, the cardinality of the actuator set that Algorithm [T] 
returns is up to a multiplicative factor of O(logn) from 
the minimum cardinality actuator sets that meet the same 
performance criterion. We note that this is the best achievable 
bound in polynomial-time for the set covering problem in the 
worst case ll24l . while © is a generalization of it (cf. ll). 


D. Approximation Algorithm for Problem (HI) 


We present an efficient approximation algorithm for Prob¬ 
lem (|I| that is based on Algorithm [T] To this end, let A be 
the actuator set returned by Algorithm [T] i.e. A G C|A| and 
(/)(A) < E. Moreover, denote as Ai, A 2 , ■ ■ ■,\n and qi, q 2 , ..., 
the eigenvalues and the corresponding orthonormal eigen¬ 
vectors of Fa, respectively. Additionally, let Am = niinjg[„] Ai 
and qM = argmax^. ig[„]r’^9i. Finally, consider a positive e 
such that ne{v^qM ^/for some c > 0. Then, 


0(A) > ^^(rA-f e/) 


Ai -f e 




> w^FA^t^ — 


neiv'^qMf 


XI. 


> — cE, 


( 11 ) 

( 12 ) 

(13) 


where we derived (El) from (fTTl ) using the fact that for any 
X > 0, l/(l-|-x) > 1 — X, while the rest follow from 
the definition of Am and qM, as well as the assumption 


nt{v^qM)"^/X^ < cE. Moreover, it is also true that 0(A) < 
E by the definition of A, and therefore from (fTSl ) we get 

< il + c)E. (14) 

Hence, we refer to c as approximation error. 

On the other hand. Am and qM are not in general known 
in advance. Hence, we need to search for a sufficiently small 
value of e so that (fl4l ) holds. One way to achieve this, since 
e is lower and upper bounded by 0 and 1 /E, respectively, is 
to perform a binary search. We implement this procedure in 
Algorithmic] where we denote as [Algorithm [1] (Fi, e) the set 
that Algorithm [T] returns, for given E and e. 


Algorithm 2 Approximation Algorithm for the Problem (U). 

Input: Upper bound E, approximation error c, bisection’s 

accuracy level a, matrices Fi, r2,..., r„, vector v. 

Output: Actuator set A. 

I 0, u -y- 1/E, e <— (I + u)/2 

while u - I > a do 

A ■(— [Algorithm [1] {E, e) 
if — v'^IFa + el)~^v > cE then 

u ■(— e 
else 
I e 
end if 

e i — (( “f ti)/2 

end while 

if v^Ff/^v — v^{Fa + el)~^v > cE then 
u i — e, e i — (/ “t“ zi)/2 

end if 

A -f- [Algorithm [1] {E, e) 


Note that in the worst case, when we first enter the while 
loop, the i f condition is not satisfied and as a result, e is set to 
a lower value. This process continues until the if condition 
is satisfied for the first time, from which point and on, the 
algorithm converges, up to the accuracy level a, to the largest 
value e of e such that u^F^^u — u^(rA -f el)~^v < cE\ 
specifically, je — ej < a/2, due to the mechanics of the 
bisection. Then, Algorithm |2] exits the while loop and the 
last if statement ensures that e is set below e so that 
— t;^(rA -f el)~^v < cE. The efficiency of this 
algorithm for Problem dH) is summarized below. 

Theorem 2 (Approximation Efficiency of Algorithm |2] for 
Problem dJ])). Denote as I* the cardinality of a solution to 
Problem © and as A the selected set by Algorithm^ Then, 

^ S ^|A|) 

v'^F-^\ <{l + c)E, (15) 

< r, (16) 

F = 0(logn + log^ + log-^-f^). (17) 

We remark that as e 0, 0(-) —>■ u^(r(.) -f eE)~^v. 
















Fig. 1. A 5-node integrator chain. 


Therefore, for any solution A° to Problem (HJi and e small 
enough, |A°| > Z*; to see this, note that as e —>■ 0 for any A°, 
(/)(A°) ^ u^(rA°+e/)“^r' < < E, since also e > 0, 

i.e. any A° is a candidate solution to Problem 0, which 
implies |A°| > I*. Therefore, for e small enough we have 
|A|/|A°| < |A|/Z*. Hence, ( fT^ is written as |A|/|A°| < F, 
that is, the worst case bound (O holds also with respect 
to the cardinality of any solution to dTJi; and as a result, the 
best-approximation properties of Algorithm [T] are inherited by 
Algorithm |2] as well. 

IV. Examples and Discussion 


bound and with the best possible actuator set for the given 
transfer, which is {1,4}, as verified in the following 

= 1.5425 • = 5.8675 • 10^ 

= 401.7997, = 6.2889, 

= 2.7445 • 10^ 

Moreover, note that although node 1 is selected as an actuator, 
in this case its corresponding input signal is zero. Thus, one 
may choose not to implement an actuator at this node, at the 
expense, however, of losing the overall network controllabil¬ 
ity. This observation motivates the analysis of dTJl when no 
controllability constraint is placed on the end actuator set. 

Finally, by setting E large enough in Algorithm |2] so that 
any actuator set respects this energy bound, we observe that 
only node 1 is selected, as expected for the satisfaction of the 
controllability constraint. 


We test the performance of the proposed algorithm over 
various systems, starting with an integrator chain in Sub¬ 
section IIV-AI and following up with Erdos-Renyi random 
networks in Subsection IIV-BI 


A. The Case of an Integrator Chain 


We first illustrate the mechanics and efficiency of Algo¬ 
rithm |2] using the integrator chain in Fig. [T] where we let 


A = 


-1 0 0 0 0 
1-10 0 0 
0 1-10 0 

0 0 1-10 

0001-1 


We first run Algorithm |2] with E <— u^E^j^gjU and a,c 
.001 and examine the transfer from a;(0) ■<— [0,0,0, 0,0]^ to 
x(l) ■<— [1,1,1,1,1]^. The algorithm returned the actuator 
set {1,4}. As expected, node 1 is chosen, and this remains 
true for any other value of a;(l), since for a chain network to 
be controllable, it is necessary and sufficient that node 1 be 
actuated. Additionally, {1,4} is the exact best actuator set for 
achieving this transfer. This is true because using MATLAB® 
we can compute 


= 5.2486 • 10®,u'^r}'^i2}^ = 2-0860 • 10^, 

= 159.9369, = 159.1712, 

= 2.1086 ■ 10^ 

Hence, node 1 alone does not satisfy the upper bound E, 
while not only satisfies this bound, but it also takes 

the smallest value among all the actuators sets of cardinality 
two that induce controllability. Therefore, {1,4} is the best 
minimal actuator set to achieve the given transfer. 

Next, we set x{l) ■<— [0,0, 0,1,0]^ in Algorithm |2] which 
led again to the selection {1,4}, as one would expect for any 
transfer that involves only the movement of the fourth node, 
while controllability is desired. In other words, even though 
we chose E •«— Algorithmic respected this energy 


B. Erdos-Renyi Random Networks 

Erdos-Renyi random graphs are commonly used to model 
real-world networked systems ll25l . According to this model, 
each edge is included in the generated graph with some 
probability p, independently of every other edge. 

We implemented this model for varying network sizes n, as 
shown in Fig. |2] where the directed edge probabilities were 
set to p — 21og(n)/n, following lH]. In particular, we first 
generated the binary adjacency matrices for each network size 
so that every edge is present independently with probability p, 
and then we replaced every non-zero entry with an independent 
standard normal variable to generate a randomly weighted 
graph. 

To avoid the computational difficulties associated with 
the integral equation Q we worked with the controllability 
Gramian instead, which for a stable system can be efficiently 
calculated from the Lyapunov equation AG-\-GA^ = —BB^ 
and is given in closed-form by 

pOO 

G= 

Jto 

Using the controllability Gramian in (|4|i corresponds to the 
minimum state transfer energy with no time constraints. There¬ 
fore, we stabilized each random instances of A by subtracting 
1.1 times the real part of their right-most eigenvalue and 
then we used the MATLAB® function gram to compute the 
corresponding controllability Gramians. 

Next, we set xq to be the zero vector and Xi the vector of 
all ones. We also set c ■<— 0.1 and a ^ 1. Finally, for each 
instance of n we first computed the corresponding lower bound 
of E so that (Uli is feasible, v^Gfy^v, and then run Algorithm|2] 
for E equal to kv^Gy^v, where k ranged from 2 to 2^^. 

The number of selected actuator nodes by Algorithm |2] for 
each n with respect to k is shown in Fig. |2] We observe that as 
k increases the number of actuators decreases, as one would 
expect when the energy bound of (ll]l is relaxed. In addition, we 
notice that for k large enough, so that (ll]l becomes equivalent 
to the minimal controllability problem of IT], the number of 









Fig. 2. Number of selected actuators by Algorithm [2| in Erdos-Renyi networks of several sizes n and for varying energy bounds E. For each n, the values 
of E are chosen so that the feasibility constraint m of Problem is satisfied: Specifically, for each set of values for n and k, Algorithm [2 was executed 
for E ■<— kv^G{n)^^v, where G(n)v is the controllability Gramian corresponding to the generated network of size n with B set to be the identity matrix. 


chosen actuators is one, as it was generally observed in m 
for a similar set of simulations. 

V. Concluding Remarks 

We introduced the problem of minimal actuator placement 
in a linear system so that a bound on the minimum control 
effort for a given state transfer is satisfied while controllability 
is ensured. This problem was shown to be NP-hard and to have 
a supermodular structure. Moreover, an efficient algorithm was 
provided for its solution. Finally, the efficiency of this algo¬ 
rithm was illustrated over large Erdos-Renyi random networks. 
Our future work is focused on investigating the case where no 
controllability constraint is placed on the end actuator set, as 
well as, on exploring the effects that the network topology has 
on this selection. 


Appendix A 

Prooes oe the Main Results 
A. Proposition [7] 

Note that since ||u ||2 = 1 and ui, U 2 , ■ ■ •, f'n-i are an or¬ 
thonormal basis for the null space of v, for any unit vector q of 
dimension equal to that of v it is (v^qY' (lY = 1- 

Next, assume that A ^ C|a| and let k be the corresponding 
number of non-zero eigenvalues of Fa. Therefore, k < n—1. 
Moreover, denote as Ai,A 2 ,...,A„ and qi,q 2 ,... ,qn the 
eigenvalues and orthonormal eigenvectors of Fa. We get 




Xj + e 




A,+e2 


n — k 


> V. 


Since e < l/u and A: < n — 1 we have a contradiction. ■ 
B. Proposition \2\ 

We first prove that u^(rA + el)~^v is supermodular. With 
similar steps one can show that vf (T/s. + for any 

i € [n — 1], also is. Then, the proof is complete, since the 
class of supermodular functions is closed under non-negative 
linear combinations. 

Recall that ^^^(rA + el)~^v is supermodular if and only 
if —n^(rA + el)~^v is submodular, and that a function h : 
V t-A R is submodular if and only if for any a € V the function 
ha ■ V \ {a} '-A R, where ha{A) = h{A U {a}) — h{A), is a 
non-increasing set function. In other words, if and only if for 
any Ai C A 2 C V \ {a} it holds true that ha{Ai) > ha (A 2 ). 

In our case, /ia(A) = -v'^{TALi{a} + e-f) ^ + v'^{Ta + 
el)~^v. Therefore, take any Ai C A 2 C V \ {a} and denote 
accordingly = A 2 \ Ai. Then, we aim to prove 

- '*^^(rAiU{a} + elGv + U^(rAi -f el)~^v > 

~ r!^(rAiU'E>U{a} + +'y^(rAiUX> + 

To this end and for 2 : e [0,1], set f(z) = ^^^(rAi + zVt> + 
Fa -f el)~^v, and g{z) = u^(rAi -f zTxi + el)~^v. After 
some manipulations the above inequality can be written as 
/(l)-/(0)>5(l)-5(0). 

To prove this one, it suffices to prove that df /dz > dg/dz, 
\/z £ (0,1). Denote Li{z) = Fai + zTx> -f F^ -f el and 
L2{z) = Fai -f^Fx)-Pel. Then, the df /dz > dg/dz becomes 

Li{z)~^Tt,Li{z)~^v < L2{z)~^Tt'L2{z)~^v, ( 18 ) 


e 
























where we used the fact that for any A y 0, B 0, z € {0, 1), 
£{A + zB)-'^ = -{A + zB)-^B{A + zB)-\ 

To show that this holds, hrst observe that 
both Li{z) and ^ 2 ( 2 ) are full rank. Thus, 
Li{z)~^) = p{t\ 1 '^L 2 {z)~^) = p{t]1'^) and, as a 
result, = 7^(^^/"L2(z)-l) = n{T]/'^) l20l. 

Hence, if u ^ 7?.(r^^), then (fTsT i holds trivially. Otherwise, 
if u S 7^(rp^), then 3'0 such that v = T^'^v and (fTST l is 
written equivalently 

v'^t]/^L2{z)-^T-dL2{z)-^t]!,'^v. (19) 

To prove ( fT9] l, it is sufficient to show that V 2 G [0,1] 

v)1'^L2{z)-^Tt,L2{z)~^^v^- ( 20 ) 

To this end, hrst observe that Li{z) ^ L 2 {z). This implies 
L 2 {z)~^ h Li{z)~^ II 20 I and, as a result. 

Now, since for any 0 ^ A ^ B, ^ B"^ ll20l . the previous 
inequality gives (1201 ). ■ 


C. Theorem Q] 

We hrst prove ((3), (0 and (fTOl i. and then ([7]). First, let 
Aq, Ai,... be the sequence of sets selected by Algorithm [T] 
and let I be the smallest index such that <p{Ai) < E. Then, 
A; is the set that Algorithm [1] returns, and this proves ([8]). 

Moreover, from ll22l . since for any A G V, h(A) = 
—(j){A) + 4){%) is a non-negative, non-decreasing, submodular 
function (cf. Proposition |2l, it is guaranteed for Algorithm [T] 
that 


I 

I* 


< 1 + log 
= 1 + log 


h{V) - /i(0) 
/i(V)-MAi-i) 
ne~^ — (j){V) 


Now, I is the hrst time that (j){Ai) < E, so (f>{Ai_i) > E. 
This implies (0. Moreover, observe that 0 < (j){V) so that 
from ([^ we get f < 1 -f log[ne“^/(£' — (/)(V))], which in 
turn implies (fl^ . On the other hand, since 0 < e < 1/E and 
<j){Ai) < E, Proposition [T] is in effect, i.e. (|7]) holds true. ■ 


D. Theorem \2} 

The hrst and the third statements follow directly from 
Theorem [T| For (fTsT l. hrst note that when Algorithm |2] ex¬ 
its the while loop and after the following if statement, 
u^FT^u — u^(rA -f el)~^v < cE. Additionally, u^(rA + 
el)~~v < (/)(A) < E; and as a result, (fTSl) is implied. Finally, 
for (fTTl i. note that u^F^^u — u^(rA + el)~^v < cE holds 
true when e is of the same order as n{v^ qm)^/( cX^E). Then, 
loge“^ = 0(loen -I- log l/fcAmFl), since < 1, 

which proves (1171 ) through (TlOl l. ■ 
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